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Abstract. Between 1996 and 2003 we have obtained 226 high resolution spectra of 16 stars in the field of the young 
open cluster NGC 6913, to the aim of constraining its main properties and study its internal kinematics. Twelve 
of the program stars turned out to be members, one of them probably unbound. Nine are binaries (one eclipsing 
and another double lined) and for seven of them the observations allowed to derive the orbital elements. All but 
two of the nine discovered binaries are cluster members. In spite of the young age (a few Myr), the cluster already 
shows signs that could be interpreted as evidence of dynamical relaxation and mass segregation. However, they 
may be also the result of an unconventional formation scenario. The dynamical (virial) mass as estimated from 
the radial velocity dispersion is larger than the cluster luminous mass, which may be explained by a combination 
of the optically thick interstellar cloud that occults part of the cluster, the unbound state or undetected very wide 
binary orbit of some of the members that inflate the velocity dispersion and a high inclination for the axis of a 
possible cluster angular momentum. All discovered binaries are hard enough to survive average close encounters 
within the cluster and do not yet show sign of relaxation of the orbital elements to values typical of field binaries. 

Key words. Binaries: spectroscopic - Stars: early type - ISM: bubbles - Open clusters and associations: general 
- Open clusters and associations: individual (NGC 6913) 



1. Introduction 

This is the second paper of a series devoted to the results 
of a long term, high resolution spectroscopic study of early 
type members of young open clusters, trapezium systems 
and OB associations. The aims of this series are discussed 
in Paper I (Munari and Tomasella 1999). 

NGC 6913, the topic of this paper, is a young open 
cluster harboring O-type members and lying close to the 
plane of the Galaxy (a = 20' 1 23 m 9, S = +38°32' (J2000); 
I = 76° 92, b — +0°61). Despite appearing in the Messier 
catalog as M29, few papers in literature deal with it, 
furthermore showing some disagreement in the results. 
Cluster distance is reported to be 2.2 kpc by Morgan 
and Harris (Q.956) and Massey et al. (1995), 1.5 kpc by 
Joshi et al. (|P3g3jl . and 1.1 kpc by Hoag et al. $TWl\ . 
while Tifft [flUSgjl suggested that NGC 6913 is indeed the 
results of two separate groups of stars, one at 1.6 kpc 
and the other somewhere between 1.9 and 2.4 kpc. The 
mean and differential reddening span a range of values too: 
< E B -v >=0.78, according to Joshi et al. 

(JEHU), < E B -v >=0.71 and A£ B _y=1.82 for Wang and 
Hu 1|2UUU|1 . and < E B -v >=1.03 following Massey et al. 
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(1995). Similarly, estimated ages span from 0.3—1.75 Myr 
of Joshi et al. ((T^jl to 10 Myr of Lynga ijTWJl . 

The internal and galactic kinematics of NGC 6913 has 
not been so far investigated in literature. The cluster ra- 
dial velocity used by Hron (1987) in modeling the rotation 
curve of the Galaxy, —25 km sec -1 , was assembled by 
scanty literature data that apparently missed all bright- 
est cluster members, and is largely off our much more 
accurate and representative — 16.9(±0.6) km sec -1 value 
(see sect. 3.2). Internal kinematics and binary content of 
NGC 6913 are unknown because no detailed radial ve- 
locity study of its members has been ever pursued, and 
proper motions investigations (Sanders 119731 Dias et al. 
2002J are not deep and accurate enough for a firm mem- 
bership segregation over a wide range of magnitudes, do 
not cover all candidate members and do not allow resolu- 
tion of the internal kinematics. 

In this paper we aim to look in more details to 
NGC 6913 general properties (like astrometric member- 
ship, photometry, reddening, distance, mass and age) and 
to present and discuss the results of our extensive spec- 
troscopic study of NGC 6913 based on 226 high resolution 
spectra monitoring of 16 stars in the field of the cluster 
over the time span 1996-2003. These observations are used 
to constrain the internal velocity dispersion, the cluster 
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Fig. 1. Finding chart for NGC 6913 program stars. 



galactic motion, the individual rotational velocities, and 
the internal kincmatical and evolutionary status of the 
cluster. Spectroscopic orbits are calculated for the discov- 
ered binary stars. 

2. Spectroscopic observations 

Table 1 summarizes the main properties of the 16 se- 
lected program stars, and Figure 1 provides a finding 
chart for them. The program stars have been spectro- 
scopically observed over the period 1996-2003 with the 
1.82 cm telescope and Echelle+CCD spectrograph of the 
Astronomical Observatory of Padova at Asiago (Italy). 
Table 2 provides the journal of observations. The instru- 
mental set-up, spectra extraction and calibration, accura- 
cies, etc. are identical to Paper I and the reader is referred 
to it for details. 

2.1. Spectral classification and radial velocities 

Scanty information exists on the spectral classification of 
the program stars. Wang and Hu {2223 derived spectral 
types from low resolution spectra (5.3 A/pix) covering 
the range 4200-6900 A. Kazlauskas and Jasevicius (1986) 
obtained photoelectric photometry in the Vilnius system, 
that we have converted into spectral types using the red- 
dening free color parameters Q defined by Strayzis (1977) 
appropriate for the Ry = Ay / Eb-v—S.6 reddening law 
that applies to NGC 6913 according to Johnson (1962). 
We have also derived spectral classification of the pro- 
gram stars using our Echelle spectra, classified against 
the Yamashita et al. (1977) spectral atlas. Even if spec- 
tral classification of Echelle spectra has to be carried out 
with care (lines to be compared normally fall on different 
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Fig. 2. The program stars on the reddening corrected V , 
(B — V)o diagram using Eb-v from Table 3 and V, B — V 
photometry from Table 1 for Ry — Av/Eb-v =3.6 ap- 
propriate for NGC 6913 according to Johnson (1962). The 
isochrone for solar metallicity and 5 Myr is from Bertelli 
et al. (1994) and it is scaled to m — M=10.5 (dotted line), 
m — Af=11.0 (solid line) and m — M=11.7 (dashed line). 
It is evident how the cluster distance cannot be well con- 
strained. In this paper we adopt a 1.6 kpc distance. 

Echelle orders), nevertheless the resulting spectral types 
look quite reasonable, and, given the far superior spec- 
tral resolution and high S/N, also possibly more accurate 
than those of Wang and Hu (2000 ) . The three estimates 
of the spectral type are compared in Table 3. The last two 
columns of the table give the reddening and distance when 
Fitzgerald (1970) intrinsic colors and our spectral classifi- 
cation are compared to V, B — V photometry in Table 1. 
The positions of the program stars on the reddening cor- 
rected HR diagram are shown in Figure 2. 

Radial velocities from individual observations (here- 
after referred to as epoch radial velocities) of the program 
stars are given in Table 4. For O and B type program 
stars they rest on individual measurement of Hel and Hell 
lines. For the other, cooler program stars the radial veloc- 
ities come from measurement of the metallic absorptions 
lines (mainly Fe I, Mg I, Ti II). The radial velocities of the 
Be program star ^5 pertain to the emission lines, which 
completely fill helium and hydrogen absorption lines. 

2.2. Binaries and orbital solutions 

About half of the program stars have turned out to be 
spectroscopic binaries. Table 5 summarizes the bary cen- 
tric velocity, the membership and the binary status based 
on epoch radial velocities in Table 4. Table 6 gives the 
spectroscopic orbits computed for all the binary stars but 
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Table 1. Program stars. The first four columns give our identification number (cf. finding chart in Figure 1), and 
that assigned by Hoag et al. (1961), Sanders (1973) and Kazlauskas and Jasevicius (1986). V and B — V are Tycho-2 
Vt and (B — V)t transformed into Johnson system following Bessell (2000) prescriptions. U — B is the median of the 
measurements by Massey et al. I)1995f) . Joshi et al. (1983) and Hoag et al. (1961). Star #10 is reported as a short 
period variable by Pena et al. 2001. 
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#6 and 16, which are clearly binaries but the available 
radial velocities are not enough to determine the orbital 
period and thus to allow to derive an orbital solution. 
Therefore the RV© quoted for stars #6 and 16 in Table 5 
is the mean of the measurements, not the barycentric ve- 
locity, and the two velocities tend to differ with increasing 
eccentricity and paucity of measurements. Consequently, 
the RV© of stars #6 and 16 quoted in Table 5 which differ 
by slightly more than 3<r from the cluster mean velocity 
cannot be considered as a firm indication that stars #6 
and 16 are field stars. 

Similarly to Paper I, the spectroscopic orbits have 
been obtained with a Fortran code written by Roger F. 
Griffin (Cambridge University) and adapted to run under 
GNU/Linux by us. 

Program stars #1-7 have been observed also by Liu et 
al. l|1989irK)91[) who reported some epoch radial velocities 
for them. Such data appear affected by large errors for the 
O and B stars (program stars #2-7), which make them 
useless in our analysis. They are instead in good agree- 
ment with our velocities for star #1, much cooler having a 
spectral type F0 III. The reason for the poor quality of the 
Liu et al. radial velocities of hot stars probably lies in the 
shortness of the wavelength range they observed (~150 A) 
and by the fact that it is dominated by H<5, which we ig- 
nored in our analysis given the Balmer progression and its 
excessive scatter compared to the much more performing 
Hcl and Hell lines. It is also worth noticing that Liu et 
al. did not recognized star #7 as double lined, in spite 
having observed it at orbital phase 0.66 when the velocity 
separation between the components is ~140 km sec -1 (cf. 
Figure 3) and therefore outstanding. 



The spectroscopic orbits of stars #2 and 11 in Figure 3 
and Table 6 are to be considered quite preliminary, given 
their small amplitude, high eccentricity and limited num- 
ber of observations. Further observations are obviously en- 
couraged for these two stars. Photometric observations of 
the double lined star #7 are in progress to the aim of 
constraining the orbital inclination and derive individual 
masses, and they will reported elsewhere when completed. 

2.3. Rotation velocities 

Rotational velocities for the program stars are given in 
Table 5. They have been derived from Hel lines for stars 
#2, 3, 4, 6, 7, 8, 10, 12, 15, 16 and Fel lines for the remain- 
ing ones, following the numerical relations for the Asiago 
Echelle spectrograph calibrated in Paper I (its Figure 7). 
No rotation velocity is derived for the Be program star #5 
because all Hel lines are badly affected by emissions. The 
correspondence of the rotational velocity scale between 
Hel and Fel lines (which we have been forced to used in 
all program stars with a spectral type later than B) has 
been carefully checked on a grid of Kurucz rotationally 
broadened spectra we have calculated on purpose. 

From the spectral classification in Table 3, the stel- 
lar radii over the HR diagram as tabulated by Straizys 
and Kuriliene (1981) and the observed V ro t sini projected 
rotation velocity, we have derived in the last column of 
Table 5 the projected rotation period (P ro t /sin?) for the 
binaries with an orbital solution in Table 6. The projected 
rotation period is obviously an upper limit to the true ro- 
tation period. Compared to the orbital period in Table 6, 
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Table 2. Journal of observations. D is the dispersions 
(A/pix) at Ha (0.19 corresponding to unbinned spectra, 
0.38 to 2 x binned spectra), and A A is the wavelength cov- 
erage. The last column gives the program stars observed 
in each given run. 



date D AA (A) star # 



1996 


.06. 


.28 





.38 


4100- 


-6700 


1,2,4,5,6,7,10,12,14 


1996 


.06. 


.29 


0. 


.38 


4100- 


-6700 


3,8,9,11,13 


1996 


.08. 


.01 


0. 


.38 


4100- 


-6700 


1,2,3,4,8,9 


1996 


.08. 


.02 


0. 


.38 


4100- 


-6700 


5,6,10,11,12 


1996 


.08. 


.03 





.38 


4100- 


-6700 


1,2,3,4,5,8,9 


1996 


.08. 


.04 





.38 


4100- 


-6700 


6,10,11,12 


1996 


.08. 


.20 


0. 


.38 


4100- 


-6700 


3,6,7,10,11 


1996 


.08. 


.21 


0. 


.38 


4100- 


-6700 


1,2,4,5,8,9,12 


1996 


.08. 


.28 


0. 


.38 


4100- 


-6700 


1,2,4,5,6,12 


1996 


.08. 


.29 





.38 


4100- 


-6700 


3,7,8,9,10,11,13,14 


1996 


.09. 


.06 





.38 


4100- 


-6700 


1,2,3,4,5,6,7,8,9,12,14 


1996 


.09. 


.07 





.38 


4100- 


-6700 


5,10,13 


1997 


.07. 


.25 


0. 


.38 


4100- 


-6700 


1,5 


1997 


.07. 


.26 


0. 


.38 


4100- 


-6700 


2,3,4,6,7,8,9,10,11,12 


1997 


.07. 


.29 


0. 


.38 


4100- 


-6700 


13,14 


1997 


.07. 


.30 





.38 


4100- 


-6700 


11 


1997 


.08. 


.10 





.38 


4100- 


-6700 


2,3,4,7,8,9,10 


1997 


.08. 


.11 





.38 


4100- 


-6700 


1,5,6,11 


1997 


.08. 


.12 


0. 


.38 


4100- 


-6700 


12 


1997 


.08. 


.13 


0. 


.38 


4100- 


-6700 


1,4,7,11,12 


1997 


.08. 


.14 


0. 


.38 


4100- 


-6700 


1,4,7,11,12 


1997 


.08. 


.15 





.38 


4100- 


-6700 


11 


1998 


.09. 


.09 





.38 


4500- 


-9500 


1,4,7,12 


2000 


.09. 


.12 





.38 


4500- 


-9500 


2,3,6,8,13,14 


2000 


.09. 


.13 


0. 


.38 


4500- 


-9500 


2,3,6,8,13,14 


2000 


.09. 


.14 


0. 


.38 


4500- 


-9500 


2,3,6,8 


2000 


.09. 


.16 


0. 


.38 


4500- 


-9500 


2,3,4,5,6,13,14 


2000 


.09. 


.17 





.38 


4500- 


-9500 


1,2,6,7,14 


2001 


.07. 


.27 





.19 


4500- 


-9500 


2,3,4,7,10 


2001 


.07. 


.27 





.38 


4500- 


-9500 


6 


2001 


.07. 


.28 


0. 


.38 


4500- 


-9500 


7,10,15 


2001 


.07. 


.28 


0. 


.19 


4500- 


-9500 


2,3,4,6 


2001 


.07. 


.29 


0. 


.19 


4500- 


-9500 


2,3,4,6 


2001 


.07. 


.29 





.38 


4500- 


-9500 


7,15 


2001 


.07. 


.29 





.19 


4500- 


-9500 


2,3,4,6 


2001 


.07. 


.30 





.19 


4500- 


-9500 


1 


2001 


.09. 


.09 


0. 


.38 


4500- 


-9500 


11,12 


2001 


.09. 


.09 


0. 


.19 


4500- 


-9500 


1 


2001 


.10. 


.29 


0. 


.19 


4500- 


-9500 


1,2,3,4 


2001 


.11. 


.02 





.19 


4500- 


-9500 


1,4,6 


2001 


.11. 


.29 





.19 


4500- 


-9500 


2,3 


2001 


.11. 


.30 





.19 


4500- 


-9500 


1,3,7,10 


2001 


.12 


.28 


0. 


.38 


4500- 


-9500 


3,7,15,16 


2002 


.01. 


.02 


0. 


.38 


4500- 


-9500 


7,15 


2002 


.01. 


.03 


0. 


.38 


4500- 


-9500 


16 


2002 


.01. 


.04 





.38 


4500- 


-9500 


3,7,16 


2002 


.01. 


.06 





.38 


4500- 


-9500 


15,16 


2003 


.05. 


.21 





.19 


4500- 


-9500 


1,5,16 


2003 


.05 


.22 





.19 


4500- 


-9500 


16 



it can be used to infer about the co-rotation status of the 
binaries. 



Table 3. Spectral types of the program stars from Wang 
& Hu ( 2000.), from photometry in the Vilnius system by 
Kazlauskas and Jasevicius (1986) transformed by us into 
spectral types following Strayzis (1977), and from classi- 
fication of our Echelle spectra against the Yamashita et 
al. (1977) reference spectral atlas. The last columns give 
the Eb-v (from Fitzgerald 1970 intrinsic colors) and the 
spectro-photometric distances for our spectral classifica- 
tion and the photometry in Table 1. 







spectral type 






# 


WH00 


Vilnius 


our 


Eb-v 


d (kpc 


1 


F0 III 


Fl V 


F0 III 


0.11 


0.2 


2 


B0 I 


Bl la 


B0 I 


0.97 


2.0 


3 


B0 I 


Bl la 


B0 I/II 


1.06 


1.4 


4 


07 II 


Bl la 


09 lb 


0.92 


2.3 


5 


B0 Ille 




B0 Ille 


1.07 


1.2 


6 


B0 II 


B2 la 


B0 V 


0.93 


1.1 


7 


B4 II 


B3 V 


B5 IV 


0.26 


1.6 


8 


B0 III 


B0 III 


B0 Ib/II 


0.99 


3.0 


9 


A4 V 


F0 V 


A4 V 


0.23 


0.3 


10 


B6 II 


B2 V 


B5 V/IV 


0.41 


1.2 


11 




A5 V 


A5 III /ll 


0.39 


1.9 


12 




B2 V 


B2 V 


0.26 


5.3 


13 


G6 III 




G5 III 


0.00 


1.1 


14 


A8 II 


A4 V 


F7 IV 


0.00 


0.7 


15 


Bl II 


B3 lb 


Bl III 


1.38 


1.6 


16 




Bl II 


B0 V 


0.83 


2.8 



Table 4. Example of the Table containing the epoch ra- 
dial velocities (and their errors) for the program stars, 
available in full in electronic form only. 



# 5 



HJD 


RV 


o 


err 


50263 


.424 


-19. 


.8 


1. 


.5 


50296. 


.514 


-19. 


.5 


1. 


.2 


50299. 


.487 


-18. 


.1 


0. 


.9 


50316 


.566 


-18. 


.8 


0. 


.9 


50324. 


.479 


-19 


.3 





.9 


50333 


.434 


-19 


.8 





.9 


50333 


.559 


-19. 


.8 


0. 


.9 


50655. 


.500 


-20. 


.5 


0. 


.9 


50671. 


.566 


-19. 


.9 


1. 


.7 


51895. 


.302 


-26 


.3 





.9 


52781. 


.402 


-18. 


.9 


1 


.0 



Star #1 is an SB1 eclipsing binary and therefore the 
sini projection factor converge toward unity, which allows 
a direct comparison between rotation and orbital periods, 
the former being twice longer. The lack of synchronicity 
could be related to the primary evolving away from the 
main sequence and the time scales of the two processes. 



C. Boeche et al.: Kinematics and binaries in young stellar aggregates. II. NGC 6913=M29 



5 



krn sec~ 



1 , , , , 1 , II 1,1111,1 II 1 1 




#1 


— — 1 — 1 — 1 — 1 — — 1 — 1 — 1 — h 


1 ' ' H ' ' ' 1 _ 

V #2 \ ~ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




#3 


1 1 1 1 1 1 1 1 h 

- -/* 


— | — I — I — I — I — | — I — I — I — I — | — 




\ / 


^-| 1 1 1 1 1 1 1 h 

: o 


— 1 — 1 — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — 1— 

o : 


: r*. cf °H 

\ / 


• A 9ti 

A> Jr O TO 

/ \/ 1 " 


- T °/ 


M / \ M " 

Y OP 
& #7 I 


: — i — i — i — i — i — — i — i — i — h 


— | — I — I — I — I — I — I — j — - 


- i , , , , i , , , , 


1 - 




•n #12 : 


i \ •/ 


\ •/ - 




i , , , . i , , , , i 



0.0 



0.5 1.0 1.5 

orbital phase 



2.0 



Fig. 3. Orbital solutions for the binary program stars (cf. 
Table 6). 



Given the masses estimated from the spectral type and 
the amplitude of radial velocity variation, star #7 prob- 
ably has a high inclination too, possibly being eclipsing 
itself. The rotational velocity in Table 5 pertains to the 
B5 IV primary, the measurement of the secondary being 



Table 5. Heliocentric radial velocity (with its standard er- 
ror) of the program stars, binary status and cluster mem- 
bership according to radial velocities, and projected ro- 
tational velocity (with its standard error). The radial ve- 
locity of the binary stars is the barycentric velocity from 
the orbital solutions in Table 6. The last column gives the 
projected rotation period P ro t /sin i of the solved binaries 
to be compared with the orbital period. 



# 


RV 


binary 


member 


V ro t sin i 


P r0 i/sini 




(km s" 1 ) 






(km s" 1 ) 


(days) 


1 


4-8.6±0.5 


yes 


no 


26±1 


5.0 


2 


-17±6 


yes 


yes 


128±3 


10 


3 


-18.5±0.6 


yes 


yes 


84±2 


15 


4 


-15.0±0.6 


yes 


yes 


103±2 


10 


5 


-19.4±0.3 


no 


yes 






6 


-25.6±0.9 


yes 


(yes) 


253±5 




7 


-18±3 


yes 


yes 


46±5 


3.8 


8 


-15.5±0.4 


no 


yes 


67±2 




9 


-5.6 ±0.5 


no 


no 


22±2 




10 


-14.8±0.7 


no 


yes 


64±7 




11 


0±3 


yes 


no 


41±6 


11 


12 


-17±6 


yes 


yes 


165±7 


1.5 


13 


-14.6±0.2 


no 


yes 


0±2 




14 


-19.9±0.4 


no 


yes 


11±3 




15 


-16±2 


no 


yes 


115±4 




16 


-27±4 


yes 


(yes) 


325±20 





too uncertain given the difference in brightness. The ro- 
tation and orbital periods are quite close, and in view of 
the uncertainties at play, the primary in star #7 looks 
synchronized. 

Stars #3 and 12 are evidently not co-rotating, be- 
cause the projected rotation period is at least several times 
shorter than the orbital period in Table 6, and working on 
sini can only enlarge the difference. For the remaining bi- 
nary stars #2, 4 and 11, no conclusion can be drawn about 
the co-rotation status, the projected rotation period being 
longer than the orbital one. 



2.4. Cluster membership 



Sanders 1)19731 hereafter S73) has published an astromet- 
ric investigation of 228 stars in the field of NGC 6913, 
identifying 92 possible members. He has however used only 
one plate pair, with an epoch separation of just 22 yr, 
with moreover the first epoch plate "severely blackened 
by the moon" . Consequently, noting the too large frac- 
tion of detected members among the measured stars, he 
warned that the cluster separation from the field is not 
satisfactory, and that the member/non-member status he 
assigned may be frequently in error. S73 limiting mag- 
nitude is V=13.8, with a completness limit not fainter 
than 1^=13.0 that corresponds to 1.2 M Q on the main 
sequence of NGC 6913. Dias et al. $M)% hereafter D02) 
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Fig. 4. Comparison between Dias et al. (2002) and 
Sanders (1973) membership data for the 24 stars in com- 
mon. 



have used proper motions from the Tycho-2 catalog to ac- 
complish the astrometric member segregation, following 
the analytical approach of Sanders (1971). They have in 
common 24 stars with S73. As Figure 4 shows, for 6 of the 
24 common stars the membership status of S73 and D02 
are in disagreement, and for the remaining 18 there is a 
fair agreement. 

A firm membership segregation is required for any 
kincmatical investigation of the cluster, and to achieve the 
best possible result, the astrometric data should be com- 
plemented by radial velocities and placing of the program 
stars on the HR diagram. 

The radial velocity distribution of the program stars 
is presented in Figure 5, where the cluster grouping 
at —16.9 km sec -1 is evident, with a dispersion of 
cr=2.0 km sec -1 . 

Table 7 summarizes the membership status according 
to the astrometric investigations of S73 and D02, the pho- 
tometry presented and discussed by Crawford et al. (1977) 
and Joshi et al. (1983), the spectral classification of Wang 
and Hu (2000), our reddening free HR diagram of Figure 2 
and the radial velocities in Table 5. Their combination 
provide our final, adopted membership reported in the 
last column of Table 7. 



3. Internal kinematic of NGC 6913 

3.1. A foreground cloud hiding part of the cluster 

As reviewed in the Introduction, all previous investiga- 
tions of NGC 6913 agree on the large differential and to- 




Fig. 5. Radial velocity distribution of the program stars. 
A Gaussian with center at RV = —16.9 km sec -1 and 
<7=2.0 km/sec fits the member stars. Shaded elements rep- 
resent stars which radial velocity differs by more than 3cr 
from the cluster mean velocity. 



charts it is evident how NGC 6913 lies close to a very thick 
interstellar cloud that seem to hide part of it. 

Our and literature estimates about the cluster distance 
converge toward a 1.6 kpc value. Studies of the interstellar 
extinction toward NGC 6913 (Crawford et al. 1977, Neckle 
and Klare 1980) agree on a steep increase of the extinction 
at 1 kpc, about halfway the distance to the cluster, as if a 
single, major cloud is responsible for the majority of the 
extinction all the way to NGC 6913. 

In Figure 6 we compare the same field centered on 
NGC 6913 as seen on Palomar POSS-II blue charts and 
by IRAS satellite at 100 /im. The north-east quadrant is 
clearly deprived of stars in the optical image, while it is 
bright in the far infrared, a clear sign of thick dust absorb- 
ing in the optical and emitting in the IR. Star counts from 
USNO-A2.0 (stars detected in both blue and red POSS-I 
prints) and near-infrared 2MASS survey (stars detected 
both in J as well as H and K bands) strongly support 
the argument of a strong foreground interstellar extinc- 
tion crossing the field of NGC 6913 and increasing steeply 
toward the north-east quadrant. 

We expect this foreground optically thick cloud to hide 
from view part of the cluster, even if the cluster center 
seems confidently identifiable with the grouping of the 
massive O and B stars seen in the optical. 

3.2. Cluster mass 

There are no published estimates of the NGC 6913 total 
mass. Two static approaches are considered in this sub- 
section, a dynamical one is investigated in sect. 3.4. 

A lower limit to the cluster mass is obtained by adding 
the mass appropriate to the spectral type of known mem- 
bers (luminous mass). Taking Wang and Hu (2000) spec- 



tal reddening affecting the cluster. Inspecting the Palomar tral types of all S73 probable members classified by them 
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Table 6. Orbital solution for the binary stars we discovered in NGC 6913. For binaries #6 and 16 the epoch radial 
velocities in Table 4 did not allowed a determination of the orbital period and therefore the derivation of the orbit. 
The errors are given in parenthesis in units of the last digit. The last raw gives the r.m.s. deviation of the solution 
from the observed radial velocities. 



orb. element program star # 







1 


2 


3 


4 


7 


11 


12 


p 


(days) 


2.70466(1) 


1.7075(1) 


697(7) 


3.51042(5) 


3.4588(1) 


1.56842(4) 


4.0350(5) 


e 




0.0 


0.8(7) 


0.21(5) 


0.0 


0.35(2) 


0.4(2) 


0.13(9) 


K 


(km sec -1 ) 


60.1(7) 


26(6) 


13.5(6) 


50.9(8) 


87(2), 125(5) 


16(5) 


79(6) 


7 


(km sec" 1 ) 


+8.6(5) 


-17(6) 


-18.5(6) 


-15.0(6) 


-18(3) 


0(3) 


-17(6) 


To 


(+2450000) 


1132.626(4) 


1362.2(1) 


652(25) 


892.40(1) 


923.96(6) 


782.3(1) 


668.7(3) 


UJ 


(deg) 


0.0 


187(14) 


134(12) 


0.0 


282(7) 


283(49) 


344(26) 


asini 


(10 6 km) 


2.24(2) 


0.4(17) 


126(6) 


2.45(4) 


3.84(9) 


0.3(1) 


4.3(3) 


f(m) 


(km sec -1 ) 


0.061(2) 


0.001(9) 


0.16(2) 


0.048(2) 


0.19(1) 


0.0005(5) 


0.20(4) 


r.m.s. 


1.41 


1.09 


1.36 


1.72 


1.57 


1.49 


1.57 



and calibration into masses from Straizys and Kurilicne 
(1981), it results 

M cl ~ 700 M (1) 

S73 warned that a sizeable fraction of his members could 
be spurious. We assume here that they compensate for 
those fainter than the S73 completness magnitude, and 
therefore 700 M Q is taken as a fair indication of the total 
luminous mass of the cluster. 

Assuming that the members of NGC 6913 distribute 
according to the Salpeter (1955) law N(m) = C to -2 - 35 
offers another possibility to estimate the cluster mass. It 
seems fair to assume that all O and B type cluster mem- 
bers have been detected and recognized as such in the 
Wang and Hu (2000) spectral survey of NGC 6913. They 
are 31 in total, spanning the range between 6 and 67 M . 
This allows to estimate the constant C 
mi 

31 = C I m _2 - 35 dm — > C = 489 (2) 
J& 

which provides an initial total stellar mass for the cluster 
amounting to 

/■120 

M cl = 489 / N{m) m dm ~ 3100 M Q (3) 

J 0.08 

distributed in about 10 4 member stars. This is an upper 
limit to current cluster mass, because (a) stars became un- 
bound early in the cluster evolution when it started to lose 
gaseous mass blowed away by the energetic winds of the 
first massive stars that formed, (6) the relaxation mecha- 
nism leads to evaporation of the lighter members, and (c) 
the mass function appears to flatten toward lower masses 
(cf. Briceho et al. 2002). The Salpeter's and observed 
mass functions are compared in Figure 7. In sect. 3. 5 and 
Figure 8 will later show how there are possible evidences 
that the cluster is at least partially relaxed. Supposing 
that the associated evaporation of members affects stars 
fainter than S73 completness magnitude, it is found that 

/■67 

M cl = 489 / N{m) m dm ~ 1000 M Q (4) 

J 1.2 



distributed in 350 members, comparable to the luminous 
mass in Eq.(l). 

3.3. Tidal and half-mass radii 

Cluster member venturing on orbits wider than the cluster 
tidal radius have a fair chance to become unbound due to 
the action of the gravitational field of the Galaxy. The 
cluster tidal radius is defined as (cf. Binney and Tremaine 
1987): 




where Rq, M g are the radius of the cluster galactic or- 
bit and the galactic mass contained within (i?G=8.3 kpc 
and M G =9.5 10 10 Mq). Inserting M c ;=700 M Q provides 
Rt ~ 10 pc, corresponding to 20' at the estimated 1.6 kpc 
cluster distance. On photographs, the cluster tidal radius 
is conventionally taken as the distance from the cluster 
center at which the stellar density drop to field value, es- 
timated by Lynga, (1987) to be - 10' for NGC 6913, which 
corresponds to 4.7 pc at the cluster adopted 1.6 kpc dis- 
tance. For R t =4.7 pc Eq.(5) would provide a much lighter 
cluster mass, M c /=50 M Q . Eq.(5) describes a statistical 
relation of equilibrium between cluster and galactic poten- 
tials that can be realized only on times scales long enough 
for perturbation fluctuations to average. Given the very 
young age of NGC 6913 (corresponding to just 2% of its 
galactic orbital period), the conditions supporting Eq.(5) 
can hardly be matched, and a M c i estimated via Eq.(5) 
cannot be trusted. For the same reasons, the R t estimated 
from star counts has to be preferred to the value estimated 
using Eq.(5). 

The half-mass radius (Rhm) is a useful quantity fre- 
quently used in N-body simulations. We have estimated it 
by measuring the radius that contain half of the 700 Mq 
luminous mass of Eq.(l), which turned out to be 8', corre- 
sponding to R; lm =3.7 pc at distance of 1.6 kpc. Following 
Binney and Tremaine (1987) the virial radius can be ex- 
pressed as 
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Table 7. Cluster membership of program stars from liter- 
ature according to astrometric, photometric, and spectral 
type criteria, and ours based on spectrophotometric paral- 
laxes and radial velocities. The last column gives our final 
membership status obtained by merging the results of the 
various criteria. S73 = Sanders 1)1973[1 . D02 = Dias at al. 
(2002:), C77 = Crawford et al. (1977), J83 = Joshi et al. 
(1983), WOO = Wang and Hu (2000). 





astrom. 


photom. 


spctr 


ours 




# 


S73 D02 


C77 J83 


WOO 


dist 


rv 




1 

1 


y 


ii 


n 


ii 


n 


n 


n 




o 


y 


y 


y 


y 


y 


y 


y 


Y 
i 


3 


y 


y 


y 


y 


y 


y 


y 


Y 


4 


y 


y 


y 


y 


y 


y 


y 


Y 


5 


y 


y 


y 


y 


y 


y 


y 


Y 


6 


y 


y 


y 


y 


y 


y 


(y) 


Y 


7 


y 


y 


ii 


y 


y 


y 


y 


Y 


8 


y 


y 


y 


y 


y 


y 


y 


Y 


9 


ii 


ii 


y 






ii 


ii 


N 


10 


y 


y 


ii 


y 


y 


y 


y 


Y 


11 


u 


ii 


ii 






y 


ii 


N 


12 


n 


ii 


y 






ii 


y 


(N) 


13 


y 


y 


ii 


y 


y 


y 


y 


Y 


14 


y 


y 


ii 


y 


y 


ii 


y 


Y 


15 


y 


y 




n 


y 


y 


y 


Y 


16 




y 








y 


(y) 


Y 



R v 



Rh; 



0.4 



9 pc 



(6) 



3.4. Internal velocity dispersion and dynamical mass 

The mean radial velocity of NGC 6913 members is 
RVg = -16.9 ± 0.6 km/sec (7) 



with an observed dispersion of a o b s —2.0 km sec -1 . The 
average uncertainty of the radial velocity of individual 
members is <7j ns t r =1.0 km sec -1 , which leads to an in- 
trinsic radial velocity dispersion of NGC 6913 members 
amounting to 



a rv = 1.70 km sec 



(8) 



the three quantities being related as a rv = y/ o 1 obs — cr'f nstr . 
Many of the physical quantities for NGC 6913 derived in 
the following critically depend upon this number. We are 
confident that our observations, reductions and measure- 
ments are state-of-the-art and involve so many high resolu- 
tion observations of such a large number of stars over a so 
long period of time, that doing much better with 2m-class 
telescopes is hardly feasible. On the other hand, the high 
competition to access larger instruments prevents them to 
be assigned on programs like the present one that require 
such a massive telescope usage (even if equipped with 
MOS devices). Undoubtly, observations of more member 
stars over longer period of times would be useful to better 
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Fig. 7. Mass function of the cluster NGC 6913 as it results 
from the spectral classification of the member stars. The 
dashed line is a fit with a Salpeter law N(m) = 489 m" 2 35 
to O and B stars (M >6 M Q ). 

constrain the intrinsic radial velocity dispersion, but this 
docs not appear as an easy task. Therefore, we believe 
that in absence of deeper/better efforts, our measurement 
of the intrinsic radial velocity dispersion of NGC 6913 is 
worth some dynamical considerations that we develop in 
the rest of this paper. 

The virial theorem links the cluster mass M within the 
radius R to the velocity dispersion in the form 



-y GM 



2R V 



3a, 



(9) 



For a virial radius of 9 pc, the virial cluster mass becomes 

(10) 



M, 



G 



4 10 4 Mi 



There is a large difference between luminous and virial 
total mass for the cluster (a factor of sixty, 7 10 2 vs 
4 10 4 M©). Could it be intra-cluster mass ? Wang and 
Hu (2000) have derived a differential reddening across the 
cluster amounting to A£?B_y=1.82 mag. Even suppos- 
ing the effect is solely caused by intra-cluster extinction 
(from above discussion the major contribution to differ- 
ential reddening should actually come from the optically 
thick, foreground dust cloud) the amount of dust required 
to produce it is of the order of 4 10 3 M , derived assuming 
a simple spherical symmetry, constant density, a standard 
dust-to-gas ratio [N{H) = 5.810 21 i?B-v atoms cm -2 , 
Kilian 1992) and that a Ai?B_y=1.82 mag differential ex- 
tinction is accumulated as result of the line of sight cross- 
ing the whole cluster. Even under these assumptions the 
corresponding amount of intracluster material is far from 
filling the gap between the luminous (i.e. Eq.l) and dy- 
namical mass (i.e. Eq.10). 

Possible explanations to solve such a discrepancy could 
be any combination of the following effects: 
(?) the large and optically thick interstellar cloud discussed 
in sect. 3.1 hides from view a significant portion of the 
cluster, causing an underestimate of the tidal radius and 
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Fig. 6. Comparison of a 30'x30' field from POSS-II (top left) and IRAS 100 fim wavelength image (bottom left) 
centered on NGC 6913. Stellar surface density of the same stellar field is obtained from USNO-A2.0 optical catalog 
(top right) with minimum and maximum density of 0.5 and 4.7 stars/arcmin 2 (clearest and darkest zone) respectively 
and grey levels in step of 0.25 stars/arcmin 2 , and from 2MASS infrared catalog (bottom right) with minimum and 
maximum density of 6.9 and 14.2 stars/arcmin 2 (clearest and darkest zone) respectively and grey levels in step of 0.5 
stars/arcmin 2 . 



luminous mass from star counts. Deep infrared imaging at 
K band and longer wavelengths could test this scenario; 
(ii) the cluster is still relaxing and some of the stars 
considered as members are actually unbound leaving 
the cluster at velocities just larger than the escape 
one, inflating the apparent dispersion of radial veloci- 
ties. An example could be star #6 which is a member 



according to both astrometric investigations and com- 
bined photometric+spectroscopic criteria, and lies pro- 
jected close to cluster center sporting one of the ear- 
liest spectral types (B0 V). Its radial velocity is well 
determined in — 25.6(±0.9) km sec -1 , which is however 
8.7 km sec" 1 away from the cluster mean velocity of 
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However, this is actually an upper limit to the mass below 
which the members tend to leave the cluster due to evap- 
oration. In fact, if the cluster core seems already relaxing 
in spite of the very young age, this could not be yet the 
case for the lower density outer regions of the cluster, with 
the relaxation being a process moving outward from the 
cluster center on timescales longer than the NGC 6913 
age. Only devoted, deep photometric investigations can 
address what is the turning mass in NGC 6913 for which 
lighter members are already experiencing evaporation and 
what is its radial dependence. 



Fig. 8. Member stars of NGC 6913 with mass smaller 
(left) and larger (right) then 6 M . Massive stars show 
a clear clustering while the lighter seems randomly dis- 
tributed. 



— 16.9(±0.6) km sec 1 . The difference exceeds 4er b s . 
escape velocity from NGC 6913 is 



= y/l2a% v = 5.9 km sec 



The 



(11) 



supporting the possibility that star #6 is actually escap- 
ing. Kroupa et al. (2001) N-body simulations of young 
clusters confirm that the observed dispersion of radial ve- 
locities could be biased by unbound members; 
(Hi) some cluster members could be components of very 
wide binaries, which orbital motion amount to just a few 
km sec -1 and cannot be detected by observations span- 
ning less than some decades; 

(iv) a cluster angular moment larger than zero, with the 
axis pointing away from the line of sight, could pass un- 
noticed given the restricted number of observed stars and 
still contribute a fraction of the observed radial velocity 
dispersion. A cluster rotation having a projected radial 
component of — 1.5 km sec -1 at half mass radius could 
be enough to bring the virial mass in agreement with the 
luminous mass. 

3.5. Evaporation of members 

Lighter cluster members seem evaporating from NGC 6913 
as the comparison in Figure 7 between the observed and 
Salpeter mass function supports, and as mass segregation 
in Figure 8 suggests. 

In fact, for O and B type members the mean mass is 
m ~ 22 Mq and being a rv = 1.70 km sec -1 the radial 
velocity intrinsic dispersion, the mean kinetic energy of O 
and B type members is 
I 3 

E c = -m» 2 = -Mai, = 1.9 1 45 erg 
2 2 



(12) 



If the cluster is relaxing toward energy equipartition, this 
is the mean kinetic energy to which its members are aim- 
ing, and there must be a value of the mass below which 
the velocity exceeds the escape value and the stars tend 
to evaporate. Such limiting mass mu m follows from 

' " ~~ mi im = 5.5M (13) 



3.6. Crossing time, relaxing time, mass segregation 

Crossing time and relaxing time are theoretical quantities 
which play a relevant role in cluster dynamics. In fact they 
are closely related to the mass of the cluster, its dynamical 
status and the number of members. The crossing time is 
related to virial radius and cluster mass by (v from Eq.7) 



tr.r 



2R,, 



GM C , 



(14) 



and the relaxing time relates to number of members, clus- 
ter mass and half mass radius as 



H 3 



N 



GM cl 81og(0.4iV) 



(15) 



(from Spitzer 1987, derived for the case of globular clusters 
with similar mass members). 

In the case of the luminous mass, M c i —700 M Q 
and N=92 (cf setc. 2.4), it is t cr =A5 Myr and t re t ax =30 
Myr, significantly longer that the estimated cluster age 
(-5 Myr). For the Salpeter's mass M cl -3100 M and 
A=10 4 (cf. Eq. 3), the crossing time reduces to t cr =20 
Myr while the relaxation time goes up to t re i ax — 650 Myr, 
both still larger than cluster lifetime. 

NGC 6913 appears relaxed, at least in its core (where 
the relaxation time is expected to be shorter given the 
higher mass density) as mass segregation in Figure 8 
shows. It is worth to note that concentration of massive 
stars toward cluster center is observed in some clusters to 
be present since their birth and not as a the result of purely 
dynamical evolution (i.e. the Orion trapezium system, see 
Hillenbrand & Hartmann, I1998J) and that some N-body 
simulations (Portegies Zwart et al. 200^1 show mass seg- 
regation to happen in clusters over ages shorter then the 
canonical t re i ax . So, the apparent concentration of heav- 
ier members toward the center of NGC 6913 could also be 
due to mechanism other than dynamical relaxation. 

3.7. Binaries 



IS G 



-m iim < sc = £c = 1.910*°erg 



The dispersion of velocities in NGC 6913 
\fio rv =2.9 km sec -1 . Binaries with orbital velocities 
faster than it tend to survive close encounters, those or- 
biting slower risk ionization (Kroupa 2000). An orbital 
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velocity of 2.9 km sec -1 corresponds to an orbital period 
of 10 6 - 3 days for a binary with a total mass of 5.5 M , 
which would go clearly undetected in the course of a 6-yr 
long monitoring like ours. Such a binary would have an 
angular separation of 0.4 arcsec (within detection thresh- 
old of current observational techniques from the ground) 
which would rise to 4.4 arcsec for members with a total 
mass of 67 M Q and decrease to 0.07 arcsec for members 
with a total mass of 1 M© . 

All binaries detected in this investigation appear 
strongly bound, not ionizable by close encounters with 
other cluster members, and quite probably primitive (the 
short cluster age argue against a capture scenario). Their 
large eccentricities and non synchronous orbits indicates 
how far they still are from tidal circularization of the or- 
bits and locking of the rotation and orbital periods which 
characterize the field binaries. Wider, more ionizable bi- 
naries are beyond the realm of spectroscopy, and would be 
profitably searched for by high spatial resolution imaging. 

4. Desiderata 

A deep and wide field photometric investigation of 
NGC 6913 and surrounding field would be a good start- 
ing point to better constraint the total mass, tidal radius 
and drop in the luminosity function of the cluster, and to 
address the large discrepancy between observed luminous 
mass, integrated IMF mass and the observed virial mass. 
At the cluster distance and reddening, UBVRI photome- 
try complete to V = 20 will map all cluster members more 
massive than 0.8 M©, thus venturing well into the realm 
of masses that should be already evaporating from the 
cluster. Such a photometric investigation, which is highly 
encouraged, should extend over a radius of not less that 20' 
from the cluster center and should include JHKL bands 
to overcome the very strong differential extinction caused 
by the foreground interstellar cloud discussed in sect. 3.1. 
Once members lighter than those here investigated will be 
firmly identified, a study of their radial velocity distribu- 
tion would add a good deal of constraints to the kinemat- 
ical status and evolution scenario of NGC 6913. 
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Table 5. Epoch radial velocities of the program stars (full 



table given in electronic form 



only). 



